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ABSTRACT 

We present a three-dimensional simulation of the corona of an FK Com-type rapidly rotating G 
giant using a magnetohydrodynamic model that was originally developed for the solar corona in order 
to capture the more realistic, non-potential coronal structure. We drive the simulation with surface 
maps for the radial magnetic field obtained from a stellar dynamo model of the FK Com system. This 
enables us to obtain the coronal structure for different field topologies representing different periods 
of time. We find that the corona of such an FK Com-like star, including the large scale coronal 
loops, is dominated by a strong toroidal component of the magnetic field. This is a result of part of 
the field being dragged by the radial outflow, while the other part remains attached to the rapidly 
rotating stellar surface. This tangling of the magnetic field, in addition to a reduction in the radial flow 
component, leads to a flattening of the gas density profile with distance in the inner part of the corona. 
The three-dimensional simulation provides a global view of the coronal structure. Some aspects of 
the results, such as the toroidal wrapping of the magnetic field, should also be applicable to coronae 
on fast rotators in general, which our study shows can be considerably different from the well-studied 
and well-observed solar corona. Studying the global structure of such coronae should also lead to a 
better understanding of their related stellar processes, such as flares and coronal mass ejections, and 
in particular, should lead to an improved understanding of mass and angular momentum loss from 
such systems. 

Subject headings: stars: coronae - stars: activity - stars: magnetic field 



1. INTRODUCTION 

A massive object that rotates very quickly and whose 
magnetic field flip-flops can significantly affect the sur- 
rounding environment. Such an object is the single, 
late-type giant G star FK Comae Berenices (FK Com; 
HD 117555), which is the epony mous prototype for th e 
so-called "FK Com-type stars" ( |Bopp fc Stencel||1981[ ). 
The FK Com stars are observed to have extremely fast 
rotation accompanied by enhanced stellar activity and 
are thought to have evolved from coalesced binaries. 

FK Com itself has been observed since the 1960s, and 



at the beginning of the 199Qs Jetsu et al.| ( |1991 [19931 
presented a long-term study or the magnetic activity of 
the system. They found that over a period of 25 years 
the location of the most concentrated spot activity has 
been moving back and forth between the same two lon- 
gitudes, which are separated by 180 degrees. They nick- 
named this activity behavior a "Flip-Flop" of the stellar 
magnetic field. 

In the past two decades, many observations of the 
FK Com system have been undertaken to better resolve 
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Korho nen et al. 

pbooj ) have performed spectroscopic and 



loto- 



metric studies of t he evolution of surface spot activity 



on FK Com, while Korhonen & Elstner (2005) studied 
this evolution using a stellar dynamo model. ^Observa- 
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tions across th e electromagn etic spectru m, inc luding ra- 
dio bands (e.g [Hughes fc Mc Lean 19 87} [Rucin ski 1991 J 
reveal broad Ha emission (e.g Ramsey et al.'l981 ; Walter 
i2 Basri 1982a; Holtzman & Nations 1984; Kjurkchicva ^ 
[Marchcv 2005), strong UV emission in which transition 
region lines are broadened to full-width half-maxima in 
excess of 500 km s^^, which is at least t wice the pro- 



jected surface equatorial rotation velocity (Ayres et al 



2006 1 , and X-ray activity greatly elevated compared with 
normal stars at this age and stag e of evolution (e.g Wal- 
ter|1981[ [Drake et al. 2008 ; Buzasi et al |[2003| 
The rapid rotation of I' K Com is probably t 



ipid rotation of I' K Com is probably the source 
of the observed enhanced activity, while the longitudi- 
nal flip-flop of the active spot area is related to stellar 
dynamo action. The flip-flop does not necessarily corre- 
spond to a complete reversal of the magnetic field, but 
rather a reversal in the location of particular, more ac- 
tive longitudes due to a non-axisymmetric stellar dynamo 
(e.g Berdyugina 2004 ). Evidence for magnetic flip -flop 
behavior on the iSun (Berdyugina fc Usoskin 2003) sug- 
gests that this phenornenon might be a universal feature 
of any realistic, non-axisymmetric stellar dynamo. 

The magnetic flip-flop by itself is intriguing. Here, 
however, we ask what would be the coronal structure 
of such fast rotating, flip-flopping star? In the case of 
slowly rotating stars like the Sun, the structure of the 
stellar corona and the stellar wind is dominated by the 
topology of the stellar magnetic field. A common ap- 
proximation for describing both solar and stellar coronae 
is the so-called "potential field approximation" . Under 
this approximation, there is no forcing on the magnetic 
field (i.e., magnetostatic field) so the magnetic field can 
be described as a gradient of a scalar potential, and the 
three dimensional distribution of the magnetic field can 
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be obtained by solving Laplace's equation for this scalar 
potential. 

This solution requires two boundary conditions for the 
spatial domain. The inner boundary condition can be 
obtained using available surface maps of the stellar mag- 
netic field. The outer boundary condition, however, can- 
not be obtained from observations. It is set in a more 
arbitrary manner, by the requirement that at this bound- 
ary the magnetic field is purely radial. This assumption 
is reasonable if we consider that the stellar wind over- 
comes the magnetic pressure of the stellar magnetic field, 
so the latter becomes fully open above the Alfven point 
(at which the Alfvenic Mach number. Ma = u/ua = 1, 
with ua — B/ -y47r/9 being the Alfven speed) . It is conve- 
nient to define a surface ("source" or Alfvenic surface), 
which is the manifold defined by the Alfvenic points. 
This is commonly approximated using a spherical surface 
with a certain height above the stellar surface, although 
the correctness of such an approxima tion has been un- 
der debate even for the solar case (e.g [Riley et~al][2006 



Gilbert et al.|2007 ). An even more serious problem arises 



when the very basic assumption, that the forces in the 
system are negligible compared to the magnetic force, 
breaks down. The rotation period of FK Com is 2.4 



days dBopp fc Rucin ski' 198?, 'Walter fc B asri|1982b||Ko 
rhonen~erarpUnn| ' Ayrcs ct al. 2006 ) , and the equatorial 
azimuthal speed, — rS2^, in the low corona (1 — 3.5-R*) 
is about 250 — 700 km s^^ . At these speeds, the dynamic 
pressure is Pdyn = pu^ ~ 0.01 — 1 dyne cm~^. Even if we 
consider a strong field on the surface with a field strength 
of Bq = 250 G, the magnetic pressure at r = 3.5R^, will 
be of the order of Pm — Bq/{3.5^ ■ Stt) 1 dyne cm~^. 
Therefore, in these rapidly rotating systems, the ratio of 
Pdyn/Pm IS uot very small, and the potential field approx- 
imation is no longer valid. A more physical method to 
describe the coronal structure is required. 

In this paper, we describe a magnetohydrodynamic 
(MHD) simulation of the corona of an FK Com-like star 
based on the flip-flop dynamo model of t he su rface mag- 



netic field of l Elstner fc Korhonen ( 2005[) and Korhonen 
fc Elstner | ( |2005l ). We adapt a global MHD model, orig- 
3lopea 



inally developed for the solar corona, to the FK Com 
model system. This model enables us to study the global 
structure of the corona and to obtain the steady-state, 
non-potential solution that includes the stellar wind and 
the effects of rapid rotation. Another advantage of the 
MHD solution over the potential field approximation is 
that it provides the solution for the complete set of phys- 
ical parameters describing the coronal plasma: the mag- 
netic field, gas density, temperature, velocity and electric 
currents throughout the three-dimensional volume under 
consideration, while the potential field only provides in- 
formation about the magnetic field distribution. 

The paper is structured as follows. The numerical 
model and the observational constraints used in the 
simulation are described in Section [21 The results 
are presented in Section [Sj and the main findings are 
discussed in Section |4| We summarize this work in 
Section [5] 



2. NUMERICAL SIMULATION 

2.1. MODEL DESCRIPTION AND 
OBSERVATIONAL CONSTRAINTS 

For the simulation of an FK Com-like star, we us e 
the solar corona model by Cohen et al. (2007 2008), 



which is part o f the Space Weat her Modeling ij'rame 



work (SWMF) ( Toth et al. 2005 ) and is based on th e 
generic MHD BA'i'S-R- US model ( [Powell et all[l999| . 
The model is driven by surface magnetic field maps that 
are used to both determine the initial (potential) mag- 
netic field distribution as well as to scale the boundary 
conditions on the stellar photosphere. In this model, a 
stellar wind solution is generated to be consistent with 
the distribution of the surface magnetic field. 

Energy deposition in our simulations is assumed to be 
controlled by the expansion factor, fg, of flux tubes, as 
in the case of the Sun. The basic assumption is that the 
energization of each parcel of the stellar wind depends on 
the amount of expa nsion of the flux tube fr om which this 
parcel comes from. Wang & Sheeley ( 1990 1 defined an fg 
as the ratio between the magnetic flux at the source sur- 
face (where the fleld becomes radial) and the magnetic 
flux at the photosphere for a particular flux tube. Based 
on this deflnition and by fltting solar wind data, they 
derived an empirical relation for Ugw j^, where Ugw is 
the terminal wind velocity, u{r — >■ oo), for each flux tube. 
This inverse relation states that slow wind comes from 
regions of large expansion of the flux tube (the bound- 
ary between open and closed field lines), while fast wind 
comes from regions with small expansion (nearly parallel 
open field lines) . An improved empi rical formula has 
been presented by Arge & Pizzo ( 2000 ) and is known as 
the Wang-Sheeley-Arge model (WSA). This model has 
been used to predict the solar wind distribution based 
on the solar surface magnetic field and its potential field 
extrapolation. However, it does not provide any physical 
parameters other than the wind speed and the polarity 
of the interplanetary magnetic field. 

In order to obtain the complete set of physical param- 
eters from our MHD model, we use the WSA model to 
determine how much excess energy to apply in each part 
of the simulation domain such that the input wind dis- 
tribution is recovered in the MHD solution. Assuming 
the conservation of energy along a streamline (or a flux 
tube) , we can equalize the total energy of the stellar wind 
on the stellar surface and at infinity. The total energy 
of the wind far from the star is equal to the bulk kinetic 
energy of the plasma, while the total energy on the stel- 
lar surface equals to the enthalpy minus the gravitational 
potential energy: 



70 fcbTo GM^ 
R. ' 



70-1 



(1) 



where 70 is the surface value of the polytropic index, 7, 
fcf, is the Boltzmann constant, rup is the proton mass, G is 
the gravitational constant, and i?* and are the stellar 
radius and mass, respectively. To is the surface temper- 
ature; this boundary condition is a free parameter in the 
model. In principle. To is taken to be equal to the aver- 
age temperature of the stellar corona, but we emphasis 
that this is only a parameter in the model, which does 
not attempt to describe the actual coronal temperature. 



3 



The kinetic energy due to stellar rotation is omitted as- 
suming that this integral describes the radial acceleration 
of the wind and that the azimuthal component does not 
contribute to this acceleration. The Bernoulli equation 
equalizes the two ends of the streamline and in principle, 
the azimuthal component of the kinetic energy could be 
estimated and be added to the equation. However, due 
to the lack of stellar wind observations to constrain our 
model, we choose to adopt the original, solar empirical 
input. We expect the dominant effect of the fast rotation 
to appear in the solution mostly due to the rotational 
forces applied on the coronal plasma. The energy equa- 
tion in the MHD solution does include all components of 
the velocity. 

Eq. [l] enables us to relate the value of 70 to the termi- 
nal speed, Usw^ originating from this point and is known 
from the WSA model. Observations reveal that the value 
of 7 close to the Sun is found to be close to unity, while 



it is close to 3/2 in the solar wind (Totten et al. 1995 



choose to use of the solar wind input over the use of a few 
indirec t observation and gener al scaling laws for stellar 
winds ( [Wood et al.][2004l [20051 ), which hold many other 
uncertainties, in any case, the range of magnitude and 
distribution of the wind from the empirical formula in 
the WSA model can be easily modified, but this param- 
eterization is out of the scope of this paper. In addition, 
the basic assumption here is that as long as the process 
of wind acceleration is assumed to be similar to the Sun, 
the relation Usw{^/ fs) holds. Since a complete theoret- 
ical model for wind acceleration is unavailable even for 
the Sun, we argue that our approach can provide an in- 
sight about the structure of stellar coronae, despite of 
the many assumptions it holds. 

2.2. SIMULATION SETUP 

FK Com rotates so rapidly that its absorption lines 
are strongly smeared, presenting a severe challenge to 
Zeeman-Doppler imaging techniques that have been sue 



19961. This is due to the fact that close to the Sun the cessfuUy applied to other stars (e.g Donati et al. 2009 



plasma contains some amount of "turbulent" internal en 
ergy, which is released in the process of the wind accel- 
eration, so that far from the Sun the plasma is much less 
turbulent with 7 close to 3/2. T he complete description 
of this concept can be found in Roussev et al.| ( 2003 ) . 
Based on this concept, and the relation 7o(Wsjuj, we can 
determine a volumetric heating function Ej(r, 70), which 
deposits energy in flux tubes such that smaller values of 
7o (which have larger "internal" energy) , lead to a faster 
wind speed. 

Once is specified, the model solves the set of con- 
servation laws for mass, momentum, magnetic induction, 
and energy: 



dt 



V • (pu) = 0, 
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with r = 1.5 (note that the parameter 7 that defines 
the empirical energy source term E-y is always < F) until 
convergence is achieved. As in any other numerical MHD 
model, the condition of V • B = needs to be enforced 
throughout the simulation. In this mod el, the "eight 
wave" method is used for divB cleaning (Powell et al. 
19991 [T6th|[2000| . 



The necessary inputs for the model are the surface dis- 
tribution of the radial magnetic field, the boundary value 
for the density, pq, and the temperature, Tq, as well as 
the stellar radius, i?*, mass, M^,, and rotation period, 
ri*. The stellar properties used in the si mulation are 
prese nted in Table [IT (see the summary of Drake et al. 
|2008'). The two boundary conditions (po and I'o) are 
chosen so that they are slightly higher than the typical 
values used for the solar would be expected for an 

active star. One can argue that the empirical wind speed 
calculated in the WSA model is fitted to the solar case, 
and that it should not be used for stellar applications. 
However, due to the lack of stellar wind observations, we 



Donati fc Landstreet||2009 ). Consequently, surface mag- 



netic field maps ar e not available and we use the su rface 
field predicted by [Elstner fc Korhonen[ ([^005 ) and Ko- 
rhonen fc Elstner ( |2U05[ ) based on a dynamo simulation 
of the system, where the maximum value of w 250 G is 
in agreemen t with recent magneti c field observations of 



the system (Korhonen et al. 2009). The dynamo model 



predicts the temporal evolution of the stellar field, and 
for the purpose of this paper, it provides the surface dis- 
tribution for different phases of the flip-flop cycle. We 
simulate the coronal structure of the FK Com model for 
three different phases, at which there are distinct differ- 
ences in surface magnetic topology. Figure [l] shows the 
magnetic maps used for the three cases, which we call 
"Case A" , "Case B" , and "Case C" . Case A corresponds 
to the held being concentrated in two single large spots 
of opposite polarity located in opposite latitudinal hemi- 
spheres. Case C has a more equal division of field in two 
spots of opposite polarity in each latitudinal hemisphere, 
again mirrored in polarity about the equator (so the total 
number of spots is 4). Case B is intermediate between 
the two, in which there is one dominant spot in each lat- 
itudinal hemisphere, but with a second spot containing 
significant field of opposite polarity. We do not consider 
the cases differing to these in phase by 180 degrees in the 
magnetic cycle, since these would simply repeat the mag- 
netic structure of cases A-C mirrored about the equator. 
Our main focus is to study the global effect of fast rota- 
tion on the coronal structure of these models. The maps 
used here have low resolution relative to solar maps, so 
we do not need to capture small active regions like in 
the solar case, and the simulation described here does 
not require very high resolution in order to capture the 
global coronal structure. The smallest grid size near the 
surface of the star is Aa; = 3 • lO^^i?^ « 2 x 10^° cm, and 
the total number of cells used is 3 • 10^. 

We ran the simulation in the frame of reference ro- 
tatin g with the star using the local time step algorithm 
( (Cohen et al. 20081. This enables a much faster con- 
vergence tor steady state simulations. The limitation of 
this approach is that the rotational forces become too 
large at large distances when simulating a very rapidly 
rotating system, and the simulation becomes numerically 
unstable. Here, since we are interested in the low coronal 
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structure, we limit the Cartesian simulation domain to 
extend only up to ^5R^, in each direction. 

3. RESULTS 

We first compare a potential field extrapolation with 
the non-potential fields from the corresponding station- 
ary MHD simulation in Figure [2] The potential field ex- 
trapolation (left panel) has a source surface (not shown 
in the zoomed in figure) located at 2.5i?^, and a number 
of field lines, both open and closed, are shown. Field lines 
that cross the source surface are fully open, as required 
by the boundary condition mentioned in Similarly, 
field lines corresponding to the same footpoint locations 
are also shown for the non-potential field from the MHD 
solution for Case A (right panel; the other cases show 
qualitatively similar behavior, as we show below in Fig- 
ure |3|. The field hues are color-coded based on their 
behavior in the two solutions: those that are closed in 
both are colored blue, those that are open in both are in 
yellow, and those that change their topology between the 
potential and MHD solution are shown in red. The blue, 
closed field lines are all quite low-lying in both solutions, 
though are significantly more stretched out in the MHD 
solution. The behavior of the red field lines is instead 
quite striking: these represent a class of field lines that 
are stretched and twisted when a wind flow and rotation 
are imposed on the system. 

Another important difference between the potential 
field and MHD solutions is the effect of the fast stel- 
lar rotation on the geometry of the field lines. Due to 
the (infinite) conductivity of the plasma in the MHD 
solution, the magnetic field lines are frozen in to the 
plasma, which propagates radially in the inertial frame. 
The footpoints of the field lines are attached to the rotat- 
ing stellar surface. As a result, the field is wound up and 
stretched around to for m an enhan ced, compact version 
of the Parker spiral (Parker 19581. A major difference 
here is that in the solar case, the toroidal component of 
the coronal field becomes dominant beyond the Alfven 
point. In the case of rapidly rotating stars, as exempli- 
fied by our simulations, the toroidal component is strong 
inside the Alfven point, and can feedback on the system. 
This strong toroidal component cannot be predicted by 
the potential field approximation. 

The global structure of the magnetic field and the ra- 
dial velocity fields derived from the MHD simulations are 
shown in Figure [s] for Cases A (top), B (middle), and C 
(bottom). The panels show slices in the field at dif- 
ferent orientations: one to indicate the 3D structure by 
showing the solutions in the y = and z = planes 
(left), and one to show the structure in the y = plane 
containing the rotation axis (right). The stellar surface 
is shown colored with the input radial magnetic field, as 
in Figure [21 and the 3D magnetic field lines are shown 
in white. The differences in the structure of the velocity 
fields in the three cases are more apparent in the right 
panel figures, where there is a clear correlation between 
the distribution of the radial stellar wind and the surface 
magnetic topology. In Case A, we have a single domi- 
nant fast stream in each latitudinal hemisphere associ- 
ated with the main strong spots. In case C, we obtain 
three fast streams together with a weaker fast stream, all 
associated with the equally strong surface spots. In case 
B, the magnitude of the radial flow decreases, since in 



this case the surface spots are weaker than cases A and 
C. The weak outflow component of the flow introduces a 
more turbulent solution and a strong inflow component 
in the southern hemisphere. Another notable feature is 
the fact that magnetic field lines that originate from re- 
gions close to the strong spots have a shape similar to 
the one predicted by the Parker spiral (due to the strong 
outflow component), while field lines that originate from 
regions far from the strong spots (where the outflow com- 
ponent is weaker) are more likely to be affected by local 
changes of the flow. 

The stellar wind structure is determined by the topol- 
ogy of the open/closed held lines and the expansion of the 
flux tubes. In the solar case, where the solar rotation is 
not extreme, the input speed from the WS A model is well 
reproduced by the MHD solution ( Cohen et al.|200 8|. In 
the case of FKCom however, we use the empirical niodel 
to specify the acceleration of the radial wind, while we 
use the MHD model to study the effect of fast rotation 
on the coronal structure. This part cannot be captured 
by the empirical model. The difference in the stream 
structures for the three cases can be attributed to the 
different density and magnetic field structures present in 
the stationary solutions. These are illustrated in Fig- 
ure H] that shows the number density close to the star 
in the y = plane, as well as an isosurface of constant 
density, n = 1 ■ 10^ cm~'^ (green surface). It can be seen 
that the large streamers in Case B (middle panel) are 
more stretched than case A (left panel) and case C (right 
panel) . This reduces the number of flux tubes associated 
with open fleld lines (with small expansion) and, as a 
result, the radial outflow component decreases as seen in 
Figure [3j 

4. DISCUSSION 

Based on the simulation results, we emphasize two 
aspects of the results. First, the non-potential, MHD 
steady state solution is significantly different from the 
potential field extrapolation for the FK Com-like flip-flop 
dynamo system. Second, the coronal structure of such an 
FK Com-like star is expected to be qualitatively different 
from the solar corona and coronae of more slowly rotat- 
ing Sun-like stars. The former is evident from the class of 
field lines that change their topology dramatically from 
the potential field to the MHD simulation (Figure [ij, 
and the latter can be seen in the large-scale topology of 
the coronal magnetic field (Figure ^pl- This structure 
is dominated by the tangling of the rotationally-wound, 
large-scale magnetic field, unlike the solar case which is 
dominated by coronal holes and active regions. 

The rapidly rotating plasma environment acts to in- 
hibit the radial component of the flow, and the den- 
sity decrease with radial distance is less pronounced than 
would be expected in Sun-like coronae. The strong ro- 
tational component causes a toroidal "dragging" and 
stretching of the field lines, which nevertheless remain 
closed since the azimuthal forces do not overcome the 
magnetic tension as in the radial case. This stretching, 
however, stores magnetic energy in the loops and causes 
an increase in the magnetic tension, Tg = B • VB//io, 
over the entire length of the loops. The magnetic tension 
is one contributer to the Lorentz force (the other one is 
the magnetic pressure) and it has the SI units of Pa 
(momentum). 
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Figure [5] shows maps of the magnetic tension in 
nPa on the x = (left), y — (middle), and z = 
(right) planes for cases A-C (top to bottom). The mag- 
netic tension in the maps is controlled by the interplay of 
two components. The outflow component that stretches 
the field lines radially, and the azimuthal component of 
the flow that stretches the closed loops in the azimuthal 
direction. In case A, we have two main fast radial streams 
and the signature for them can be seen closer to the star, 
especially in the y = plane. Case B is characterized by 
a weaker outflow component, so that the azimuthal com- 
ponent takes over to dominate the distribution of the 
magnetic tension. In particular, the magnetic tension 
on the z = plane in this case is stronger, indicating a 
stronger stretch applied on the equatorial loops. In ad- 
dition, the magnetic tension is weaker in the southern 
hemisphere due to the lack of strong outflow in these re- 
gions. Case C is characterized by four outflows at higher 
latitudes, so the magnetic tension near the equator is 
weaker than the other two cases. 

In Figure [6] we show the distribution of the ratio of 
Tb in the MHD solution to that in the potential field, 
for Case C on the y = and z = planes (the overall 
behavior for Cases A and B is similar). The signature 
of the stretched field lines is easily visible in the equa- 
torial plane (z = 0) plots. The magnetic tension in the 
y = plane is mostly due to the radial gradient in B 
and, therefore, it is quite similar in the potential and 
non-potential solutions. The magnetic tension is highest 
near the surface, particularly near regions of spot activ- 
ity where active regions are expected to exist. These are 
regions where the loop footpoints will be subjected to 
considera ble stress due to convective motions in th e pho- 

Parker|[l974| [Fisk|[2005l ) . Such 



monly found o n rapidly rotating dwarfs (e.g Robinson 
fc Collier Came ron 1986; Collier Cameron & Robinson 
1989a|b||ColJiCT Cameron fc Woods||1992p : 



tospherc (Leighton 1964 



sites are conducive to enhanced activity and the gener- 
ation of X-ray fiares. In the simulation presented here, 
the corona reaches a steady state so that the stretched 
loops remain unchanged. It is likely, however, that time- 
dependent processes observed on the Sun, such as foot- 
point surface motions and evolving small active regions 
exist on FK Com-type stars that can trigger stellar erup- 
tions and fiares. It is also possible that these processes, 
in particular loop footpoint motions, can lead to mag- 
netic reconnections between the stretched loops them- 
selves and to major flaring activity on a global scale from 
the star. The study of these processes would be of great 
interest for rapidly rotating coronae but cannot be ad- 
dressed in the stationary simulations presented here. 

We expect that the strong azimuthal wrapping of the 
coronal magnetic field, within the Alfven radius, will be a 
general feature of the coronae of rapidly rotating stars of 
all spectral types. We also draw attention to the closed, 
wrapped magnetic field structures evident in the left pan- 
els of Figure |3] and in Figure |4] The apex of such a 
field line lies between one to several stellar radii from the 
stellar surface. This type of structure would be a likely 
candidate for hosting the "slingshot" prominences com- 



5. SUMMARY AND CONCLUSIONS 

We have performed global MHD simulations for 
an FK Corn-like system in order to study the three- 
dimensional structure of the corona of very rapidly ro- 
tating stars. We drove the model using surface synthetic 
magnetic maps produced by a stellar "fiip-flop" dynamo 
simulation of FK Com. 

Our main findings are: 

1. For rapidly rotating systems, the potential field ex- 
trapolation is inadequate to describe the structure 
of the coronal magnetic fields; 

2. The simulations show the presence of a significant 
azimuthal component in the coronal fiow that im- 
pedes the radial component (the stellar wind out- 
flow); 

3. The fast stellar rotation combined with the highly 
conductive outflowing wind generates a strong 
toroidal component of the magnetic field in the 
form of highly stretched coronal loops within the 
Alfven radius and wrapped open field lines beyond; 

4. This stretching introduces regions of large mag- 
netic tension in the stationary solutions which may 
act as sites where magnetic reconnection can be 
preferentially triggered via surface footpoint mo- 
tions and emergence of new fiux. A disconnection 
of such stretched loops might lead to major, large- 
scale stellar flares. 

The simulation presented here provides a first glimpse 
of the structure of the corona of a rapidly rotating star. 
We expect the salient features of our results to be gen- 
erally applicable to rapidly rotating dwarfs as well as 
FK Com-type stars. Of interest for future study would 
be the distribution and magnitude of the predicted wind 
mass fiux, and the influence on this of the changing coro- 
nal topology through the dynamo cycle. In addition, 
time-dependent study of how the azimuthally stretched 
loops interact and influence magnetic reconnection and 
the evolution of coronal mass ejections would be highly 
motivated. 
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TABLE 1 
Adopted Properties of FK Com . 
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Case A Case B Case C 




-260.0 -187.6 -126.0 -62.6 -0.0 64* 1S6.0 187.6 260.0 [G] 



Fig. 1. — Input magnetic surface maps used in Cases A. B. a nd C (left to right) obtained from the stellar "flip-flop" dynamo model for 
FK Com (^Elst ner fc Korhonen|2005||Korhonen fc Elstner|2005| |. These cases cover three different regimes of the dynamo cycle: dominance 
of each heinisphere by a single spot (Aj; two spots m each hemisphere of equal magnetic field strength (C); and the intermediate case in 
which there are two significant spots but with one dominant (B). 




Fig. 2. — Potential field extrapolation (left), and non-potential MHD solutions (right) for Case A. Blue field lines are closed in both 
solutions, yellow field lines are open in both solutions, and red field lines are field lines that their topology drastically changes from the 
potential field to the MHD solution. The stellar surface is colored according to the magnitude of the perpendicular surface magnetic field. 



8 




Fig. 3. — Global views of the stationary magnetic field solutions for the three cases A (top), B (middle), and C (bottom). The 3- 
dimensional magnetic field lines are shown in white. Color contours of are displayed on the y = and z = planes (left panels). The 
right panels show a similar display, but with a side view of the y = plane for clarity. The stellar surface is also shown, colored according 
to the magnitude of the surface field as in Figure [2] 
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Fig. 4. — Density structure and configuration of the steady state solutions for Cases A (left), B (middle), and C (right) close to the star. 
Color contours are of number density, and the green surface represents an iso surface of n = 1 ■ 10^ cm"^. The 3-dimensional magnetic 
field lines arc shown in white. 




Fig. 5. — Magnetic tension on the a; = (left), y = (middle), and z = (right) for cases A-C (top to bottom). 
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